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EQUATION OF STATE OF MATTER AT 
SUPERNUCLEAR DENSITY 

by 

E. R. Harrison 
Goddard Space Flight Center 


SUMMARY 

Using a vector meson model, ZeUdovitchhas shown that 
at very high densities, the pressure tends to equal e, where 
e is the energy density. By modifying Zel f doviteh T s model 
to avoid certain unrealistic simplifications, and by assuming 
that the higher masses of the meson hierarchy progressively 
predominate in the interactions as the density increases, it 
is shown that the pressure, in fact, tends to equal e/3. 
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EQUATION OF STATE OF MATTER AT 
SUPERNUCLEAR DENSITY 


by 

E. R. Harrison* 
Goddard Space Flight Center 


We shall consider briefly the possible nature of the equation of state of matter at supernuclear 
density. This is of interest in such subjects as cosmology and gravitational collapse. By super- 
nuclear density we mean a density exceeding that of ordinary nuclear matter. Since nucleons in 
nuclei are spaced approximately lfermi apart, we have in mind proper densities of p > 10 15 gm cm -3 . 
ZeUdovitch (Reference 1) has shown on the basis of a vector meson model that as the supernuclear 
density increases in an isotropic fluid, the pressure tends to equal the energy density. Our aim is 
to show that it is unlikely that Zel'dovitch’s conclusion is true. 

The equation of state can be expressed in the form 

P = {v - i) £ (1) 


where p is the pressure and e = pc 2 is the energy density. Both p and v are matrices which can 
be diagonalized for macroscopic elements of fluid. For a photon or neutrino gas, in which the 
particles have zero rest mass, 


where q = 1, 2, or 3 is the number of spatial degrees of freedom. When the gas is isotropic, 
q = 3, and therefore according to Equation 2, v - 4/3, and p - e/3. For a gas containing a class 
of particles of energy E i = y. m. c 2 , when m. is their rest mass and y. is the ratio of energy to 
rest energy, 


y. 2 ( q + 1) - 1 

qr/ 


(3) 


♦National Academy of Sciences - National Research Council, Resident Research Associate, on leave from the Rutherfcrgl High Energy 
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if the interaction energy is ignored. By regarding the energy contribution from the various inter- 
action fields as a gas of real and virtual bosons, the effective value of v for the fluid can be eval- 
uated by integrating Equation 3 over the various distribution functions of the particles. Since 
v. < (q + l )/q for all classes of particles we have, in general, 

q + 1 

1 < v ~ q“ (4) 


for all fluids. The lower limit is for a gas of zero pressure and the upper limit is for a gas con- 
taining only particles of zero rest mass. It is also noticed that the velocity of sound is 


V 



c(v“l) 1/2 < cq -1 / 2 


(5) 


in the appropriate direction, and has its maximum value of v s = c at q = 1 , as one would expect. 
In an isotropic gas, which is of main interest, Equation 4 becomes 1 < v < 4/3 and therefore p < e/3, 
and gives v s < (c/3) 1/2 . 

In spite of the complexity of the particle states at supernuclear density, our classical picture 
indicates that 1 < v< 4/3 must still hold true in an isotropic fluid. It would also seem reasonable 
to suppose that at supernuclear density v has a value close to 4/3. Using the adiabatic relativistic 
equation 


d(eV) + pdV = 0 


( 6 ) 


for a variable volume V, we have from Equation 1 that 



Assuming v is constant in a given density range, it follows that 


6 « v“ v (8) 

and if v has a value close to 4/3, then e « v~ 4/3 , as in a relativistic ideal gas. 

The relatively simple picture outlined above has been thrown into a state of confusion by 
Zel ! dovitch T s (Reference 1) suggestion that strong interactions increase the value of v, and that in 
the limit of high density in an isotropic fluid, v - 2, and therefore p^ e and v s - c. it is argued that 
the range of variation of v is 1 < v< 2, and not 1 <^£4/3 as given by the previous argument. An- 
isotropy in such a fluid is now impossible, as this would imply that the velocity of sound could 
exceed the velocity of light. Zel T dovitch T s suggestion is based on a simple model in which station- 
ary baryons interact through a vector meson field. 
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The Proca equations (Reference 2) for a vector meson (unit-spin) field are 


* 1 d 2 A 


V 2 A - 


c 2 at 2 


- \~ 2 A + — - T = o 


V 2 $ \ — —■ - K 2 $ + 477 p - 0 

c 2 dt 2 


(9a) 

(9b) 


with the gauge condition 


- 1 

V* A + cT ~cft = 0 (9c) 

where A and $ are the potentials and j and p e are the current and charge densities. These Lorentz 
invariant equations are for "hyperphotons" (vector mesons) of rest mass p = *iAc. When p - 0, 
they reduce to Maxwell's equations for photons of zero rest mass. ZePdovitch assumes static 
conditions and the Proca equations therefore become 


V 2 $ - XT 2 = 0 . (10) 

Hence, if g 2 /tic is the coupling constant, where g is the baryon "charge", the solution is the Yukawa- 
type potential 


$ = gr” 1 e“ r/ ^ . (11) 

ZePdovitch then shows that if there are n baryons per unit volume and m is the baryon mass, the 
energy per baryon is 


E 


me 2 


+ 



e~ r/A dV 


me 2 + 2rrg 2 n\ 2 


and the energy density is therefore 


e = nine 2 + 2rrg 2 n 2 \ 2 . (12) 

If it is assumed that the range X of the interaction is independent of density, it follows from Equa- 
tion 6 that the pressure is 


p 


d 


- dV ( 6V > 


27rg 2 n 2 \ 2 . 


(13) 
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Comparing Equations 12 and 13, it is seen that in the limit of high density we have p - e, and ac- 
cording to Equation 1, v therefore approaches a value of 2. The assumption that the potential (Equa- 
tion 11) is density- independent means that the range A. of the interaction at high densities is large 
compared with the interbaryon distance. The interaction therefore tends to become long-range and 
collective rather than short-range, thus accounting for v approaching the value of 2. 

The following comments must be made regarding Zerdovitch f s treatment: At supernuclear 
densities, it is physically unrealistic to ignore the high energy of the baryons imposed by the ex- 
clusion principle. In ordinary nuclear matter, the ideal gas laws provide a good starting point for 
calculating the energy shift (References 3 and 4) and potential functions are employed to obtain 
improved approximations. When the internuclear distance is less than the nucleon Compton wave- 
length, or p > 10 17 gm cm“ 3 , the idea of stationary baryons is completely unacceptable since the 

nucleons now have energies exceeding 1 GeV and there is copious production, among other things, 
of nucleons and hyperons (including their antiparticles). A further unsatisfactory feature of the 
model is that it singles out unit- spin mesons and neglects the rest of the meson multiplets. 

Let us, however, ignore the unsatisfactory nature of Zel T dovitch’s model and accept the po- 
tential (Equation 11) as a crude phenomenological description of short-range interactions. As the 
density increases (and the available energy involved in the interactions also increases), the higher 
masses of the meson hierarchy will progressively predominate in the interactions. In effect, k 
will get smaller and the short-range character of strong interactions will be preserved. Thus, 
the potential (Equation 11) can no longer be regarded as density independent. As an example, Levy 
(Reference 5) points out that the higher order potentials behave as exp(-xrA) and involve the ex- 
change of x mesons. It is to be expected that the range of the predominant interaction is comparable 
with the interbaryon distance, and therefore k is related to density by an expression of the kind 
k - an” 1/3 , where a is of the order unity. Equations 12 and 13 now become 


e - nmc 2 + 2rrg 2 cl 2 n 4/3 


(14) 


p 


277 

3 


g 2 a 2 n 4 ^ 3 


(15) 


and these equations have the merit that they yield the physically acceptable result that v has a 
maximum value of 4/3. 

Despite the complexity of the interactions in matter at supernuclear density, it is unlikely 
that our earlier arguments concerning the maximum possible value of v are in error. Although 
ZeTdovitch’s model possesses several unacceptable features, it is nevertheless interesting to 
notice that if it is applied in a slightly more realistic manner, then the maximum possible value 
of v in an isotropic fluid is also 4/3. 


(Manuscript received May 4, 1965) 
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